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ABSTRACT 

The weak gravitational lensing of distant galaxies by large-scale structure is ex- 
pected to become a powerful probe of dark energy. By measuring the ellipticities of 
large numbers of background galaxies, the subtle gravitational distortion called "cos- 
mic shear" can be measured and used to constrain dark energy parameters. The ob- 
served galaxy ellipticities, however, are induced not by shear but by reduced shear, 
which also accounts for slight magnifications of the images. This distinction is negligi- 
ble for present weak lensing surveys, but it will become more important as we improve 
our ability to measure and understand small-angle cosmic shear modes. I calculate the 
discrepancy between shear and reduced shear in the context of power spectra and cross 
spectra, finding the difference could be as high as 10% on the smallest accessible an- 
gular scales. I estimate how this difference will bias dark energy parameters obtained 
from two weak lensing methods: weak lensing tomography and the shear ratio method 
known as offset-linear scaling. For weak lensing tomography, ignoring the effects of 
reduced shear will cause future surveys considered by the Dark Energy Task Force to 
measure dark energy parameters that are biased by amounts comparable to their error 
bars. I advocate that reduced shear be properly accounted for in such surveys, and I 
provide a semi-analytic formula for doing so. Since reduced shear cross spectra do not 
follow an offset-linear scaling relation, the shear ratio method is similarly biased but 
with smaller significance. 



Subject headings: weak gravitational lensing, cosmic shear, dark energy 
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1. Introduction 



Weak gravitational lensing of galaxies has rapidly become an essential cosmological tool. 
Cosmic shear - the large-scale pattern of gravitational distortion embedded in the observed shapes 
of galaxies - is a particularly valuable effect that contains information about the growth of mass 
fluctuations as well as the Universe's expansion history and spatial curvatu re. Cosmic shear 
was first detected by multiple independent g r oups at the turn of th e century (|Bacon et al.l 12000 : 



Kaiser et al.l 120001 : IVan Waerbeke et al.l 120001 : IWittman et al.l 120001). and it could become one of 
the premier methods for studying dark energy ( Albrecht et al. 12006 ). Subsequent surveys, such 



as the 100 Square Degree Weak Lensing Survey, have shown that cosmic shear measurements 
provide cosmological parameter constraints that are generally consistent with mor e mature obser- 



vation s such as CMB anisotropy, galaxy correlations, and galaxy cluster counting (|Benjamin et al. 



2007b . This survey was a compilation of four separate surveys which, together, provided 3.5 mil- 
lion galaxies with a median redshift of about z med ~ 0.8. Ongoing surveys and those on the horizon 
are expected to provide wider, deeper multi-color images needed to extract information about dark 
energy. For example, NASA's Joint Dark Energy Mission is a planned space-based telescope that 
would survey 4000 deg 2 , measuring shapes and photometric redshifts for over a billion galaxies 
with a median redshift of z med ~ 1.5. 

As daunting as the observational systematics are for ambitious shear surveys, there are also 
important theoretical hurdles that must be overcome so that our ability to predict the shear signal 
will keep pace with our ability to measure it. Various theoretical systematics (uncertainties or com- 
plications in the predictions) have been identified and in some cases addressed. Determining how 
baryons and non-linear gravi tational clustering affect the matter power spectrum on ~ 1 Mpc//z 



scales is a primary concern (IHuterer & Takadal 12005 



J2008I). S ignificant theoretical 

systematics also include the effects of non-gaussianity (|Cooray & HulboOla) and the spurious sig 



Rudd et al. 



nals arising from the intrinsic alignment and clustering of source galax ies (ISchneider et al.ll2002l : 



Hamana et al.l 12002c iHeymans & Heavens! 120031 : iHirata & Seljakl 120041) . Computational approxi 



mations are also being investigated mo re closely as we star t to test regimes where simple first order 
shear calculations may no lo nger hold dvale & Whitehooi Icoorav & Hubooi IShapiro & Cooravl 



20061 : iDodelson et alJl2005t bodelson & Zhangll2005l) 



Ignoring the difference between shear and reduced shear is one example of a computational 



appro ximation that will not be accurate enough for upcoming cosmic shear surveys (|Dodelson et al. 



2006b . When a galaxy is weakly lensed, the change in ellipticity of its image is proportional to the 

reduced shear : 



l-K 



(1) 



where y is lensing shear and K is lensing convergence, a measure of the matter fluctuations pro- 
jected along the line of sight. The superscripts denote a "+" or "x" shear polarization transverse 
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to the line of sight. When the distortion fields K and y are much less than unity, we can expand 
(0Q) to first order in the fields, making the approximation that g a rs Y 1 , i.e. that galaxy ellipticities 
directly measure the shear field. This approximation is widely used by theorists since it holds for 
the vast majority of galaxy images in a cosmic shear survey and since it makes calculating the 
cosmic shear power spectrum and bispectrum much more manageable. Dodelson, Shapiro, and 
White (|2006l) have found that this seemingly innocuous substitution can nevertheless change the 
cosmic shear spectrum by as much as 10% on arcminute angular scales, where lensing is enhanced 
by the non-linear matter power spectrum. Thus, ignoring reduced shear in this regime could bias a 
survey's analysis, causing it to rule out a correct set of cosmological parameters. 

This paper assesses the impact of the reduced shear approximation on dark energy parameter 
determination from cosmic shear. Dodelson et al. considered source galaxies at a single redshift, 
and hence only assessed the bias that reduced shear would have on ACDM models. Here I general- 
ize to multi-redshift weak lensing analyses that constrain dark energy models by using information 
about distances and/or the growth of massive structures. The two methods I consider are weak 
lensing tomography and a shear ratio method. Tomography is more well-known: using the CDM 
paradigm combined with some dark energy model, one calculates a prediction for cosmic shear 
observables for sources in multiple redshift bins, and then the pr edictions a re compared with ob - 
servation; model parameters can subsequently be fit to the data (|Hull2002c: lAlbrecht et al.ll2006|) . 
Shear ratio methods forgo theoretical predictions of the shear signal and attempt to extract model- 
independent distance measurements by taking ratios of the she ar signal at d ifferen t redshifts. I 
will focus on the "offset-linear scaling" shear ratio method of Zh ang et al. (12005b . which is an 
alternative to the "cross-correlation cosmography" method of Jain & Taylor! ( 2003|) . 



The structure of this paper is as follows. In section 2, 1 will review some relevant weak lensing 
formalism and discuss how the next order correction to (OQ) alters cosmic shear power spectra. In 
section 3, 1 estimate the bias in dark energy parameters that various surveys would experience by 
using the reduced shear approximation with weak lensing tomography. I review the offset-linear 
scaling method in section 4 and similarly calculate how it is biased by reduced shear. 



2. The Reduced Shear Power Spectrum 



2.1. Review of the leading order result 



We begin with a review of relevant weak lens i ng equations (for a full introdu ction see e.g. 
Blandford et^ll99l[Eai^ll992l : lDodelsonll2003l : lBartelmann & Schneider! boOlk The lensing 
convergence at a particular sky position, k(6), can be expressed as the matter density contrast, 
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S(x), projected over comoving distance, X, along the line of sight: 

/■•oo 

Ki(0) = / d X 8(e X ,x)Wi(x). (2) 
Jo 

The subscript denotes a particular redshift bin i from which source galaxies have been selected, 
and Wj(x) is me lensing kernel for sources in the ith redshift bin, defined below. Future cosmic 
shear surveys will measure galaxy redshifts photometrically, estimating each redshift by observing 
the galaxy's flux in a few broad wavelength bands. The alternative is to obtain a redshift from 
the spectrum of each galaxy; although spectroscopic redshifts are much more precise, they will 
become prohibitively expensive and time-consuming for surveys of many millions of galaxies. If 
the source galaxies are binned according to their photometric redshifts, then 



y»r a (*) h d Xs d A {Xs 



with Wo = jQ.i„Hq. Here, pi(z) is the true (spectroscopic) distribution of galaxies in the ith photo- 
metric redshift bin, and nf al is the total 2D number density of galaxies in that bin. The scale factor 
relative to today is a = (1 and d A is a function that accounts for spatial curvature in open 

or closed Universes: 

( x for K = (flat) 

d A (x)=l \K\-^ 2 wi(\K\^ 2 x) for K>0 (closed) . (4) 
{ \K\-^ z wahl\K\^ z x) for K<0 (open) 

To leading order, the two shear fields y"(B) are easily related to the convergence in Fourier 
space: 

f(/) = r*(Z)iq(0, (5) 

where I is the Fourier conjugate of 6, and I am working in the small angle limit, rather than 
decomposing the fields into spherical harmonics. The trigonometric weighting functions are 

T\l) = cos(20 z ) T 2 (l) = sin(20z) (6) 

where 0/ is the angle between / and some fixed x-axis. It is useful to consider the following linear 
combinations of the shear modes, 

Ei(l) = T a (l)f(l) (7) 
Hi) = e ab T a {l)tf{l), (8) 

where e ab is the anti-symmetric matrix, £ 12 = — e 21 = 1. To leading order, the E-mode and the 
convergence have the same fourier coefficients, £,(/) = K/(Z), while the B-mode vanishes. These 
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relations do not account for small, sec ond-order corrections to the weak lensing fields such as 



lens-l ens coupling or source clustering (ICooray & Hull2002l : ISchneider et al]|2002l : lHamana et al 



2002J). The E-mode angular power spectra and cross spectra, Q ;iJ , are defined as 

(%{l)Bj{l')) = {2x) 2 &(! + f)C w 



(9) 



where 8 2 is a 2-dimensional Dirac delta function and the angle brackets denote an ensemble aver 
age. The leading order calculation of the Q ;;/ - is well known: 



Q 



;Wi(x)Wj(x)P S (k;x) 



(10) 



d A (x) 2 

where P§ (k;x) is the 3D matter power spectrum for k = £/cIa(x) at a distance x, accounting for 
the growth of structure. The above expression uses the Limber approximation, which assumes 
that the only matter density modes 8(k) contributing to the lensing signal are those modes with k 
transverse to the line of sight. Equation (flOl) is commonly used to predict the cosmic shear signal 
in order to constrain cosmological parameters. In particular, dark energy can be constrained via its 
effects on structure growth and on the geometrical functions, Wj(x)/dA(x)- 



2.2. Perturbative Correction 

Since galaxy shape measurements are sensitive not to shear but to reduced shear, we must 
consider the difference between the shear spectra in (flOl ) and the expected reduced shear spectra. 
We can do this perturbatively by expanding (OQ) about K = and keeping one more order: 

gf(e) = 7 ?(e) + (fK) l (6) + ... (id 

ff{K l ) @{K 2 ) 0(K 3 ) 

If one then recomputes the E-mode spectra, substituting g for f in ©, one recovers the leading 
order spectrum plus a correction term: 

5 (mEjV)) = T a (l)T b (l') ((r^Ul) f}(t))+ (I - V) (12) 

Under the Limber approximation, the power- and cross- spectra corrections from reduced shear 
simplify to an easily calculable formQ : 

8Q.ij = 2j ^cos(2^-2^)^.(/,/ / ,-/-/ / ;z) (13) 



The corresponding equation in lDodelson et al.l (12006) is incorrect as printed, containing spurious factors of / and 
/'. It should match (fT~3l > here with i — j. Here I use different conventions for the T a (l). 
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where I have defined the 2-redshift convergence bispectrum, 



Bfj(l h l 2 ,h; X ) = - 



dA(x) 



iWiWWjixmW+Wjix^BsihMMx) (14) 



with no sum over i or j and e.g. k\ = / 1 / cIa (x) ■ Thus, the reduced shear correction is proportional 
to a 2D convergence bispectrum, which can be written as a projection of the 3D matter bispec- 
trum, just as the leading order shear power spectrum is a projection of the m atter power spectrum. 
Due to the form of the matter bispectrum (|Scoccimarro & Couchmanll200ll) the integrand in (fT3l 
depends only on / • and therefore SCi-jj is independent of the direction of /. Under the Limber 
approximation, I find no corrections to the EB or BB spectra to 3rd order in the weak lensing fields. 

In 2nd order perturbation theory the matter bispectrum, B§, is proportional to P|; using this 
re sult as a starting point, an accurate fitting formula for B§ was obtained from N-body simulations 
by lScoccimarro & Couchmanl (|2001r) . I use their formula to compute the reduced shear correction, 
dCi-jj. I compute the linear matter power spectrum using the fitting fo rmula of 



Eisenstein & Hu 



(119991) . and the non-linear part is computed from the "Halofit" code of ISmith et al.l (120031) which 



has als o been calibrated by numerical simulations. This is the same procedure used by LDodelson et al 



(|2006|) . who demonstrated that the semi-analytic formula for the reduced shear correction (fT3l) ade- 
quately agrees with the difference found by ray-tracing N-body simulations (for all source galaxies 
at z = 1). Encouraged by t his agreement, and because other studies hav e shown that c orrections 
of order ?c 4 are negligible (|Cooray & Hull2002l : Ishapiro & Coorayl 12006c |KrauseJ|2008|), I do not 
refine the correction any further by including higher order terms in (fTT)) . 

Figures [TJ and [2] show the size of the reduced shear correction relative to the leading order cos- 
mic shear power spectrum. Unless otherwise noted, all calculations use the reference ACDM cos- 
mological model summarized in Tabled The correction is most significant toward the non-linear 
part of the shear spectrum, reaching several percent on small scales corresponding to lensing by in- 
dividual dark matter halos. Essentially, the correction reflects the ratio of the matter bispectrum to 
the matter power spectrum, which is enhanced by non-linearity; and indeed when the computation 
is repeated using only linear matter fluctuations, the correction is negligible on all scales. Thus 
it is not surprising that the size of the correction is sensitive to a parameter like the matter power 
spectrum normalization, o%, which strongly determines the matter power spectrum on non-linear 
scales. As shown in figure [TJ the correction is enhanced by a high value of Cg, and it would be 
similarly enhanced by raising the primordial spectral tilt n s . We can see that the correction could 
reach the 10% level for I = 2x 10 4 , but beyond this scale, shape noise will dominate even the 
deepest surveys of the near future. 

The reduced shear correction also increases with redshift, as shown in Figure [2l This is ex- 
pected on all scales since the lensing kernel Wj(x) increases with redshift, and the convergence 
bispectrum has an extra factor of Wj(x) relative to the power spectrum. Note that while the cor- 
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Fig. 1. — Dependence of the reduced shear correction on o$. The size of the correction relative to 
the leading order cosmic shear power spectrum is shown as a function of multipole. From smallest 
to largest, the curves assume <7s = 0.7, 0.8, and 0.9. All curves are for source galaxies in a single 
redshift bin with 0.8 < z < 1.0. 
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rection increases steadily on large scales, it quickly approaches a maximum on small scales. This 
reflects the fact that the small-scale bispectrum to power spectrum ratio drops off dramatically in 
the more linear, high-redshift Universe. In other words, light rays encounter the most non-linear 
structure at low redshifts. Again, when these computations are repeated using only the linear mat- 
ter power spectrum, the correction at I = 10,000 rises steadily with redshift instead of flattening 
out. 



3. Implications for Dark Energy Constraints from Tomography 

3.1. Methodology 

Having calculated the difference that reduced shear makes in the expected cosmic shear 
power- and cross- spectra, we would like to know to what extent ignoring this difference will 
bias measurements of cosmological parameters. This subsection focuses on the technique of weak 
lensing tomography, in which measurements of cosmic shear spectra are simply compared to the 
predictions of some model and then used to fit the mode l's parameters. To estimate the bias in a 



set of parameters, I first compute a Fisher matrix (see e.g. lDodelsonll2003l) which approximates the 
expected inverse covariance matrix for those parameters: 

^ = lE(2 i+1 )T r (c-.g.C r >.g) 05, 

where p a is the value of the ath parameter, and the C/ are the symmetric covariance matrices of 
the shear signal. I am assuming that the various modes / are uncorrected - this is an approximation 
discussed in more detail at the end of this section. The elements of Q are the expected shear signal 
plus shape noise, Ci-jj + Njj. The shape noise matrix is diagonal and accounts for the fact that the 
shear in redshift bin i is measured from a finite number density of background galaxies nf al with 
intrinsic shape variance y} ms : 

Afo = %. (16) 
n i 

If prior constraints on the parameters are included, the total Fisher matrix is simply the sum of the 
prior matrix and the survey matrix: 

where / sky is the sky coverage of the survey, expressed as a fraction of the total sky. Diagonal 
elements of the inverse matrix, (F tot ) _1 , are estimates of the forecasted variances of the parameters, 
marginalized over all other parameters. These estimates are accurate when the signal and noise 
distributions are approximately Gaussian. 
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Fig. 2. — Dependence of the reduced shear correction on redshift and multipole. Top: same as 
figure[Ubut for different redshift bins and Og fixed at 0.85. For source galaxies in a single redshift 
bin of width Az = 0.2, the curves are for bins with maximum redshifts n = 0.4, 0.6, 1.0, 2.0, 4.0. 
The correction increases with increasing redshift. Bottom: same correction as a function of zu 
From smallest to largest, the curves fix / = 200, 2820, 10,000. 
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If our predictions of the shear signal are incorrect, then we expect there to be a bias in (p a ), 
the expected best-fit parameter set from a survey using those predictions. Now from (TT5T) and (fTTT) 
it is straightforward to show that if the true cosmology is given by the parameters p l ™ e , then to first 
order, the parameter differences are given by 



= PT - (Pa) = ^V ot g £(2Z + l)Tr (V • |^ ■ Cr 1 • 



(18) 



where the 5Q are matrices whose elements are the "prediction errors", in this case, the reduced 
shear corrections 5C/ ;;/ . The concern is then whether this bias Ap a is large compared to the 
survey's confidence intervals about (p a ), as illustrated in Figure[3] 

If the parameter biases Ap a are substantial relative to the parameter confidence intervals ob- 
tained from the Fisher matrix, then we should be concerned about ignoring the difference between 
shear and reduced shear. However, 1 -dimensional biases and error bars do not provide the full 
picture: it is possible to slightly bias multiple parameters so as to cause their combination to be 
ruled out in the full multi-dimensional parameter space. In particular, I am interested in the dark 
energy equation of state parameters wq and w a , defined by 

w(z) =w + w a (\ -a) =w + w a ——. (19) 

1 +z 

To check whether these two parameters will be biased in co mbination, I have calc ulated a "figure 



of bias" (FoB) where appropriate (a similar metric is used in lDodelson et al.ll2006Q 



To obtain the FoB, I first invert the Fisher matrix F tot to obtain a covariance matrix for the 
cosmological parameters. I then extract the wq — w a block of that covariance matrix and call this 
submatrix (F w )~ l . The FoB is defined as 

FoB= (^ApaF^App) (20) 

W / 

where the sum only runs over the parameters wq and w a . Note that for a single parameter, the 
FoB is simply Ap/o(p). More generally, the FoB indicates the level at which we should expect 
a survey to rule out the correct dark energy parameter combination due to the reduced shear bias. 
With no bias, we expect (in a statistical sense) that the wq — w a confidence region will be centered 
on the true dark energy values like the solid ellipse in Figure |3] With a bias, we expect that 
confidence region to be displaced by the vector Ap a , thereby placing the true parameters on some 
error contour. Assuming that the error contours are ellipses defined by the Fisher matrix F w , the 
FoB squared is si mply the Ax 2 of the true dark energy point relative to the center of the ellipses 



([Press et alj 1 19921 see discussion of confidence intervals). Some common confidence levels are 
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Fig. 3. — Illustration of parameter bias Apj in the wq — w a plane. The star represents the true 
parameter values, and the solid ellipse represents the la error contour forecasted for an unbiased 
survey via a Fisher matrix analysis. The square and dotted ellipse represent the best-fit parameters 
and la error contour expected from a survey that ignores reduced shear. In this example, the 
best-fit dark energy parameters expected from the survey are incorrect by several "sigma". 



given here: 



FoB 
(FoB) 2 
Error contour 



1.52 2.15 2.48 3.03 3.44 4.29 
2.30 4.61 6.17 9.21 11.8 18.4 
68.3% 90% 95.4% 99% 99.73% 99.99% 



Thus if FoB = 1 .52, it means that we expect the true wq and w a values to lie on the 68.3% error con- 
tour of a survey's biased confidence region; FoB = 3.03 means that the true dark energy parameters 
would be "ruled out" at the 99% level. 

In order to estimate the impact of the reduced shear approximation on parameter fitting, I 
compute the Fisher matri x for three types of w eak lensing surveys described in the Report of the 
Dark Energy Task Force (lAlbrecht et al.l 120061) . The types of survey instruments I consider are a 
Stage-Ill 4-m class telescope using photometric redshifts, a Stage-IV large survey telescope (LST), 
and a Stage-IV Joint Dark Energy Mission (JDEM) satellite telescope. All surveys are assumed to 
collect a population of source galaxies with a redshift distribution given by 



dnP x 
~dz~ 



<*z 2 exp[-(z/zo) 



1.5i 



(21) 



The median redshift of this distribution is z me d = zoV2, and the distribution is normalized so that 
the total angular number density of galaxies is 



dz — ; — 

o dz 



gal = ^gal 



(22) 



I assume that photometric redshift measurements for all surveys are unbiased and have a Gaussian 
distribution about the true, spectro scopic galaxy re dshifts with a scatter of o(z phot ) = 0.02(1 +z) 
(for photo-z methodology see e.g. LMa et al.1120061) . To avoid degrading parameter constraints, it 
is actually more important to have a small uncertainty in the scatter rather than a small scatter 
(IMa et al.1120061) . The galaxies are divided into five photometric redshift bins, defined so as to 
have precisely equal shape noise when a(z phot ) = 0. Several authors have found that subdividing 
the redshifts further, w hich increases the shape noise in each bin, doe s not significantly improve 
parameter constraints (Hu 1999 : Ma et al. 20061 : Takada & Jain 2004 ). These and other survey 
characteristics are summarized in Tabled] 

The cosmological model used for forecasting is summarized in Ta ble |2l it is a standard set 
of parameters consistent with WMAP 5-year data (|Dunkley et al.l l2008h . Curvature and the dark 
energy equation of state are allowed to vary, and fiducially the model is flat with a cosmological 
constant (wq = — l,w a = 0). I have intentionally chosen a somewhat high-end value of Cg that is 
pessimistic, i.e. it increases the importance of the reduced shear correction, as shown in Figure [Q 
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Table 1: Characterist ics for three cosmic shear surveys suggested by the Dark Energy Task Force 
iAlbrecht et alJbood) . 



Parameter 


4-m 


LST 


JDEM 


Description 


ysky 


0.1 


0.5 


0.1 


Sky coverage of survey 


n£ al 


15 


40 


100 


Effective number of galaxies per arcmin 2 


z med 


1 


1 


1.5 


Median redshift of galaxies 


7rms 


0.25 


0.25 


0.3 


Intrinsic galaxy shape variability 


CT(/ cal ) 


0.01 


0.01 


0.001 


Multiplicative shear calibration error 



Photo-z bin edges 4-m, LST: 0,0.622,0.875,1.13,1.47,4 
JDEM: 0,0.930,1.31,1.69,2.19,4 



Table 2: Cosmological parameter set and fiducial values used to predict constraints and biases from 
weak lensing tomography. The last group of parameters are not independent of the others. The 
Planck Prior column contains la error bars derived from the Planck Fisher matrix. The Planck 
constraints on curvature and dark energy parameters contain a degeneracy that precludes finite 
error bars. 



Parameter 


Fiducial 


Planck 


Description 




Value 


Prior 




wo 


-1 




Dark energy equation of state today 


W a 







High-z change in dark energy equation of state 


&-DE 


0.72 




Dark energy density today 


a k 







Spatial curvature 


(Dm 


0.137 


0.0012 


Total matter density today, £l m h 2 


(Ob 


0.0225 


0.00017 


Baryon density today, &bh 


n s 


0.96 


0.0062 


Primordial scalar spectral index 


InA, 


-0.249 


0.018 


Log amplitude of primordial scalar spectrum 








at k = 0.05 Mpc" 1 


h 


0.7 




Hubble constant, H = lOO/i km/s/Mpc 




0.85 




Amplitude of linear matter power spectrum 








on 8/z Mpc~ ! scales today 


% 


2.3 x 10~ 9 




Scalar amplitude, <5 2 = A s x 2.95 x 10~ 9 
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this choice will make the final results more conservative. I impose priors expected from CMB tem- 
perature and polarization measurements from the Planck satellite; this is accomplished by adding 
Planck's Fisher matrix for the parameters in Table [2] to the Fisher matrix of each weak lensing 
survey. The Planck Fisher matrix was computed by the Dark Energy Task Force and was made 
available with the interactive program, "DETFast" (|Dickl|2006h . In addition, I marginalize over the 
nuisance parameters / ; cal to account for multiplicative shear calibration errors in each redshift bin: 
Q;ij ^ fi fj • These parameters are equal to unity, with errors given in Tabled] 

The shear modes included in my analysis are 10 < / < 3000; a high-/ cutoff is necessary 
since baryons and non-linear gravitational clustering make the shear fields on sm aller scales non- 
Gaussian and difficult to predict as of yet (|Hu & White! 1200 ll ; LRudd et al.ll2008|) . Since the Q ;;; - 
are rather smooth functions of /, in practice I compute sums by sampling logarithmically in /, and 
I have checked that the sampling is fine enough for the results to have converged. Although I 
have treated the shear modes as independent, non-li near clustering does in fact correlate the matte r 
density modes 8(k), thereby correlating the shear (|Hu & White! 1200 ll ; iRimes & Hamiltonll2005|) . 
These correlations can degrade parameter constraints (depending on the parameter), and so can 
non-Gaussianity; therefore the parameter errors obtained from the Fisher matrix will be lower 
limits. I numerically compute two-sided derivatives of the observables with step sizes of 0. 1 for wq 
and w a , -025 for n s and lnA^, and 5% for the remaining parameters; again, I check that reducing 
the step size does not significantly affect the results. 
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3.2. Results and Discussion 

Forecasted cosmological parameter constraints and biases for each of the three surveys con- 
sidered are shown in Tables [3j H and \5\ In the case of the 4-m class telescope, we can see that 
the dark energy biases are comparable to their la error bars. The equation of state parameters 
are biased by more than 0.5a in all cases, and the dark energy density can reach a la bias when 
the equation of state is allowed to vary (top two rows). The remaining parameters are primarily 
constrained by Planck. Their biases are not as significant, although n s and A s are biased by almost 
0.5a for a flat, ACDM Universe. The FoB is not very sensitive to a curvature prior. The fact that 
FoB=l . 1 indicates that ignoring reduced shear will displace the wq — w a error contours by a signif- 
icant (~ la) amount: enough to create an artificial tension with other observations. Furthermore, 
the biases may become more significant if combined with other data or priors not considered here. 
Since Stage-Hi surveys will operate in the near future, theoretical systematics may be dominated 
by other issues such as predicting the matter power spectrum. Nevertheless, accounting for reduced 
shear would be the safe approach, and doing so is straightforward via (fT3l . 



The biases from neglecting reduced shear are worse for a large survey telescope that would 
measure a higher density of galaxies and cover more sky than a 4-m class telescope. As shown in 
Table HI dark energy parameters from an LST are consistently biased beyond the la level and in 
some cases exceed a 2a bias. With FoB=2.5, the true wq — w a combination would lie beyond the 
2a error contour for those parameters. Assuming a flat ACDM Universe (last row) leads to large 
biases in n s and (O/,, parameters which can be used to tilt the shear spectra to mimic the reduced 
shear correction (with oo m fixed, CO/, changes the ratio of baryonic matter to dark matter, which 
alters early structure growth in a scale-dependent way). In addition, by the time such a Stage- 
IV survey is completed, cosmic shear predictions might be accurate enough to make neglecting 
reduced shear a dominant theoretical systematic error. 

Reduced shear will perhaps be an even more important consideration for the Joint Dark En- 
ergy Mission, wh ich is expected to have the lowest shape noise of any survey while probing the 



highest redshifts ([National Aeronautics and Space Administrationll2008l) . The top part of Table [5] 



shows that when using modes / < 3000, the dark energy parameter biases can exceed la. Biases 
in a>h, n s , and A s are also significant, with the spectral parameters breaking la when dark energy 
parameters are fixed. Furthermore, by the time JDEM is completed, cosmic shear predictions may 
be mature enough to exploit even the extremely small-scale modes which can be measured by a 
space telescope. The bottom part of Table [5] assumes that all modes with / < 10 4 will be used in 
the cosmic shear analysis. Even with only Planck priors, this assumption leads to very significant 
biases, including an FoB of 5.2, which corresponds to the correct dark energy parameters mis- 
takenly falling beyond the 99.99% confidence region of the survey. It makes sense that including 
high / modes magnifies the bias because the reduced shear correction increases towards large / and 
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Table 3: Estimated biases for weak lensing tomography with a Stage-Hi survey on a 4-m class 
telescope. For each parameter, I obtain its expected bias Ap a and its l-a error <j(p a ) from the 
Fisher matrix, assuming the survey measures shear multipoles / < 3000. Also shown is a dark 
energy Figure of Bias (FoB) as defined in the text. In addition to Planck priors (P), I include the 
following information where indicated: fixing Q. k = 0, fixing w(z) = wq, and assuming that dark 
energy is a cosmological constant (A). 



Priors 


FoB 






Wq 


W a 


&-DE 


a k 


(0,n 




n s 


InA, 


P 


1.1 


Ap a 


X 10 2 


-21 


46 


2.2 


-0.068 


0.0059 


0.00065 


0.053 


0.14 








) x 10 2 


24 


64 


2.1 


0.30 


0.12 


0.017 


0.61 


1.7 


p,a k 


1.1 


Ap a 


x 10 2 


-21 


48 


2.2 




0.029 


0.0015 


-0.023 


0.30 






<y(p a 


) x 10 2 


24 


64 


2.1 




0.060 


0.014 


0.52 


1.5 






Ap« 


x 10 2 


-3.6 




0.92 




0.026 


0.00021 


0.079 


0.25 






<y(pa 


) x 10 2 


4.6 




1.2 




0.060 


0.014 


0.50 


1.5 






Ap« 


x 10 2 






0.025 




0.015 


0.0031 


0.21 


0.64 






o(p a 


) x 10 2 






0.31 




0.059 


0.014 


0.47 


1.4 



Table 4: Same as Table [3]but for a Stage-IV Large Survey Telescope. 



Priors 


FoB 




Wo 


W a 


&DE 


a k 






n s 


InA, 


P 


2.5 


Ap a X 10 2 


-20 


43 


2.2 


-0.078 


-0.0047 


0.0053 


0.31 


0.26 






o(j>a) X 10 2 


8.8 


24 


0.86 


0.22 


0.10 


0.016 


0.56 


1.4 




2.5 


Ap a X 10 2 


-21 


47 


2.2 




0.025 


0.0018 


0.17 


0.32 






o(j>a) X 10 2 


8.4 


21 


.85 




0.057 


0.013 


0.41 


1.4 


v,a k , w 




Ap a X 10 2 


-2.5 




0.67 




0.038 


0.0085 


0.65 


-0.30 






o(j>a) X 10 2 


2.3 




0.51 




0.057 


0.013 


0.36 


1.4 


PQ^A 




Ap a X 10 2 






0.16 




0.0045 


0.015 


0.90 


0.56 






o(j>a) X 10 2 






0.22 




0.047 


0.012 


0.27 


1.1 
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Table 5: Same as Table[3]but for a space telescope survey like the Joint Dark Energy Mission. The 
lower rows add cosmic shear information from all modes with / < 10 4 . 



Priors 


FoB 






w 


W a 




Q k 


co m 


(Ob 


n s 


InA, 


P 


1.4 




X 10 2 


-9.4 


15 


1.5 


-0.013 


0.011 


0.0047 


0.31 


0.52 






<?(Pa 


) x 10 2 


8.7 


21 


1.0 


0.24 


0.11 


0.017 


0.57 


1.5 


PMk 


1.6 


Ap« 


x 10 2 


-9.5 


15 


1.5 




0.016 


0.0042 


0.29 


0.54 






<y(p a 


) x 10 2 


8.5 


19 


1.0 




0.054 


0.013 


0.41 


1.4 


P^,w 






x 10 2 


-2.9 




0.77 




0.025 


0.0050 


0.37 


0.29 






<?(Pa 


) x 10 2 


2.2 




0.48 




0.053 


0.013 


0.40 


1.4 






&Pa 


x 10 2 






0.20 




-0.013 


0.011 


0.61 


1.6 






<y(p a 


) x 10 2 






0.21 




0.045 


0.012 


0.36 


.97 


p 


5.2 


Ap« 


x 10 2 


-20 


31 


3.2 


-0.30 


-0.12 


0.026 


1.3 


0.072 


l< 10 4 




<y(p a 


) x 10 2 


6.9 


19 


0.72 


0.20 


0.099 


0.015 


0.50 


1.4 


p^ 


5.6 


&Pa 


x 10 2 


-27 


52 


3.8 




0.010 


0.012 


0.75 


0.31 


l< 10 4 




o{p a 


) x 10 2 


5.1 


12 


0.62 




0.053 


0.012 


0.36 


1.4 


P^,w 




Ap« 


x 10 2 


-6.7 




1.6 




0.11 


0.017 


1.5 


-2.9 


l< 10 4 




(?(Pa 


) x 10 2 


1.9 




0.38 




0.047 


0.012 


0.32 


1.2 


P^,A 




Ap« 


x 10 2 






0.41 




-0.017 


0.039 


2.4 


0.32 


l< 10 4 






) x 10 2 






0.16 




0.039 


0.011 


0.18 


.77 
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because the errors in the Cj-jj decrease as (21 + l) -1 / 2 for uncorrelated, Gaussian modes. We must 
keep in mind, though, that the shear modes on non-linear scales are increasingly non-Gaussian 
and correlated, and therefore the parameter errors calculated here are lower limits (FoBs are upper 
limits). 



4. Implications for Dark Energy Constraints Using Offset-Linear Scaling 

4.1. Background and Methodology 

As previously mentioned, cosmic shear is sensitive to large-scale structure formation and to 
the Universe's expansion history and curvature, a.k.a. the Universe's "geometry." It will be ex- 
tremely useful to be able to isolate these effects, that is, to use cosmic shear to probe only structure 
formation or only geometry. Doing so will allow for consistency checks on cosmic shear data. 
Isolating geometry from structure formation will also be useful for testing General Relativity since 
matter will cluster differently in different gravitational theories, even when the expansion histories 
match. Furthermore, predicting the non-linear clustering of matter is a difficult ongoing problem, 
which is why we must currently impose small-scale cutoffs in cosmic shear analyses; it would 
therefore be practical and conservative to extract the geometric information from a shear survey in 
a way that is not biased by our ignorance of the details of non-linear gravitational clustering. Shear 
ratio methods do just that - they use cosmic shear to measure the distance/redshift relation, and 
there fore constrain dark energy, without our having to predict the non-linear matter power spec- 
trum Jjain & TavloJ2003l : lzhang et all2005l : lTavlor et alJl2007h . Although shear ratio methods are 



useful for the reasons listed above, weak lensing tomography will generally yield tighter parameter 
constraints since it aggressively incorporates more information. 

The shear ratio method I will focus on in this paper is the "offset-linear scaling" method 



of IZhang et al.l (|2005l) . This method exploits the fact that a cosmic shear cross-spectrum can be 



written in the following form: 



B H 

Ci-ij ~ A /;/ + -^f f or Zi<Zj (23) 



with 



_ j dx aM»L (24, 

Xj J n) X 

Mi = W °J d *-^J a{x) 2 % , P ^X) (25) 



- 19- 



% = - W °l d X—^J a{x) 2 x , XPs(k;x). (26) 



Here, I am assuming a spatially flat Universe for simplicity. Equation (1231) follows directly from 
(flOl) and is not an approximation so long as galaxies are binned by their spectroscopic redshifts. 
When photometric redshifts are used, (1231) is valid so long as the galaxies in bins i and j have very 
little redshift overlap. The key thing to notice is that by (f24]> . the effective distance xf^ depends 
only on the distance to the jth redshift bin with some weighting given by the galaxy distribution 
(which is measured). All of the information about the matter power spectrum P§ is contained in the 
parameters A/ ;i and So it is easy to see that if we choose some common foreground redshift 
bin, i = f, then the following combination of cosmic shear spectra, 



(i. f\ — Ci;fa Q\fb _ (Zq ff ) 1 (Xif) 1 ni , 

= c l;fc -c hfd (xf)- l -(xfr 



has no dependence on P§. It also has no dependence on / or /, but depends only on distances 
to the background redshift bins. In other words, we have a cosmic shear observable that is only 
sensitive to geometric information. To constrain dark energy, we only have to predict the ^J ff , and 
we can (in principle) exploit even the non-linear high-/ information in a cosmic shear survey since 



,eff 
'3 

galaxy-shear cross-correlations; however, I will restrict my discussion to cosmic shear. 



Xj has no dependence on scale. It is possible to construct observables similar to r a ^(/;/) using 



The form of (|23l implies that any two cross-spectra with a common foreground redshift bin are 
different only by the scaling of (#? ff ) _1 , which is an offset-linear scaling. As Zhang et al. point 
out, the scaling relation is invalid if the redshifts from separate bins significantly overlap - this 
leads to extra terms in (1231) that do not fit the scaling relation. Subsequently, the ratio T a bcd{l',f) 
picks up a dependence on P§, thus ruining it as a purely geometric observable. The extra terms can 
be made negligible by reducing the scatter in photo-z measurements (a technical challenge) or by 
simply choosing well-separated redshift bins. Similarly, since we can only measure reduced shear 
and not shear, (|23T) should contain a term for the reduced shear correction, and again this disrupts 
the scaling relation. Since the reduced shear correction always contains a squared factor of the 
background lensing kernel, Wj(x) 2 , it necessarily produces a term proportional to (xf f )~ 2 '- 

C l;ij + 8C l;ij « A' l;i + + for Zi < Zj . (28 ) 

Hence a reduced shear cross-spectrum cannot be expressed as an offset-linear scaling relation. 

We can see the breakdown of offset-linear scaling in Figure HI which illustrates the fractional 
change in T a bcd{l\f) as a function of multipole for a few sample redshift bin combinations. With- 
out the reduced shear correction, these curves would all be equal to zero, indicating no dependence 



10. 50. 100. 500. 1000. 5000. 10000. 

1 



Fig. 4. — Fractional change in T a bcd{l\f) versus multipole, calculated from (l27l) but including 
the reduced shear correction. Curves assume JDEM parameters, and sources have been binned by 
photometric redshift with bin widths of Az = 0.25. The maximum redshifts of the background bins 
{a,b,c,d) are (3.5, 3.25, 1.25, 1) for the solid curve, (3.5, 3.25, 3.5, 1) for the long-dashed curve, 
(1.5, 1.25, 1.25, 1) for the short-dashed curve, and (3.5, 1.25, 3.5, 1) for the dotted curve. All 
curves have foreground bin / with 0.5 < z < 0.75. The dependence on / means that the reduced 
shear cross-spectra do not maintain the offset-linear scaling relation. 

on / (I have explicitly verified this as a consistency check). With the correction, the scale indepen- 
dence is clearly broken, and the solid curve shows a combination of redshift bins with one of the 
largest deviations. Thus, when the difference between shear and reduced shear spectra are substan- 
tial, T a bcd{l\f) is somewhat sensitive to the matter power spectrum and does not provide a purely 
geometric probe of dark energy. As previously mentioned, offset-linear scaling also breaks down 
when photometric redshift bins have substantial overlap, but in that case, we can remove scale 
dependence e.g. by not using adjacent bins. It is less clear how to remove the scale dependence 
introduced by reduced shear, and I propose one possibility in the next section. 

To estimate how dark energy parameters are biased by ignoring reduced shear in the offset- 
linear scaling method, I repeat the Fisher matrix analysis of section[3l The formalism is essentially 
the same, but now to make use of the offset-linear scaling of (|23l) . I include the A/.,- and 5/., as 
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nuisance parameters and marginalize over them so that the observables Q ;(J will only be sensitive 
to xf f - Clearly, there are now a large number of parameters to keep track of and equally large 
Fisher and covariance matrices that must be inverted; however, the matrices are sparse, and Zhang 
et al. describe how to manage them compactly^ . The parameters co m , 0) b , n s , and A s can now 
be dropped from the Fisher matrix analysis since they do not affect the #J ff , although I select 
the same fiducial values in Table [2] to compute the covariance matrices C/. The shear calibration 
nuisance parameters / ; cal can also be discarded since they are completely degenerate with the new 
nuisance parameters, A/ ; , and I marginalize the Planck Fisher matrix over all but the 3 dark 
energy parameters (fixing curvature), and I add the resulting matrix as a dark energy prior. This 
prior is degenerate, providing no constraints on its own; however it does improve constraints when 
combined with cosmic shear. 

Using the offset-linear scaling method, I repeat the Fisher matrix analysis for the Joint Dark 
Energy Mission. All of the cosmological parameters and survey characteristics remain the same as 
in section|3l except that instead of 5 redshift bins, I now use 16 evenly spaced bins from < z < 4. 
The purpose of the extra bins is to at least partially make up for the weaker constraints expected 
from the shear ratio method. As mentioned previously, finer binning is not proportionally helpful 
because the shape noise in each bin will be proportionally larger. For a higher photo-z scatter, 
wider redshift bins will be needed to avoid substantial spectroscopic overlap as mentioned above. 
Since no high-/ cutoff is needed, I include all modes with / < 20, 000. 



4.2. Results and Discussion 

As expected, comparing Table [6] with Table [5] shows that even with finer redshift binning 
and more small-scale information, dark energy parameter constraints from offset-linear scaling are 
weaker than those from weak lensing tomography in all cases. In fact, we cannot put interesting 
constraints on more than 2 parameters at a time with the offset-linear scaling method and cosmic 
shear alone. Although the parameter biases from ignoring reduced shear can be larger than the 
biases in Tabled they are not as significant. With w a and curvature fixed, wo and Q.de are biased by 
0.26c and 0.27c respectively when using modes / < 20, 000; these are not completely insignificant 
biases but neither are they a major worry. Throwing away some high-/ modes, where reduced 
shear is most important, does help to reduce the significance of the bias. Without modes / > 
3000, constraints on one or two dark energy parameters do not degrade much, yet the biases are 



2 For N bins, the Fisher matrix singular if we include the nuisance parameters in the highest two redshift bins (A/ ^, 
Bi:N, Ai-jst-i, and since no combination of C/jj can measure these. I remove the singularity by holding these 

parameters fixed; then Ci-nn-i must be ignored to avoid falsely adding information at these redshifts. 
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Table 6: Same as Table [3] but using the shear ratio method with a space telescope survey like the 
Joint Dark Energy Mission. 



Priors 


FoB 




Wq 


W a 






0.35 


Ap a x 10 2 


-47 


120 


4.2 


I < 3000 




<7(p ) x 10 2 


140 


360 


12 






A/? a X 10 2 


-1.8 




0.49 


l < 3000 




o( Pa ) x 10 2 


16 




4.6 






A^« x 10 2 






0.0014 


I < 3000 




0(p a ) X 10 2 






0.66 




0.52 


Apa X 10 2 


-53 


150 


3.5 


l< 20,000 




O(pa) X 10 2 


130 


340 


11 






Apa X 10 2 


3.9 




-1.1 


/< 20,000 




a(pa) x 10 2 


15 




4.1 


pa*, a 




Ap a X 10 2 






-0.019 


l < 20, 000 




<T(Pa) X 10 2 






0.66 



considerably smaller. Hence, decreasing the shape noise should increase bias by providing more 
weight to information on small scales. Increasing the sky coverage simply scales up the measurable 
number of modes in a survey, providing smaller error bars on dark energy parameters; since Ap a 
is nearly independent of / sky in the shear ratio method, Ap a /o(a p ) scales roughly as a/ / sky (not 
a perfect proportion when including the Planck prior). Thus, I find that if we expand the JDEM 
survey to / sky = 0.5 and fix w a , the 0.26a bias on wo becomes 0.59c 

A potential advantage of the offset-linear scaling method is its ability to use high-/ modes 
without our having to predict shear spectra for those modes. The importance of reduced shear 
means that in order to exploit these modes without bias, the method will have to be adjusted. 
When using weak lensing tomography, reduced shear can be accounted for by simply including 
the reduced shear correction 5C/ ;! -y, given by (fT3l . in our predictions. Unfortunately, since shear 
ratio methods do not require predictions of shear spectra, only the distances ^J ff , there is no pre- 
diction that can be corrected. To be clear, I have calculated cosmic shear spectra in this section for 
the purposes of forecasting, but to do this, I had to explicitly assume an underlying matter power 
spectrum. In practice, a shear ratio analysis will not depend on such information. One could con- 
ceivably guess the correction using a reasonable bispectrum prediction; however, for consistency, 
we should de-bias the shear ratio method in a way that does not involve calculating the reduced 
shear correction from theory. One option is to add an extra set of nuisance parameters, such as 
D\. i in d23l . to account for the deviation from an offset-linear scaling relation. However, this will 
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likely degrade dark energy constraints to the point of being uninteresting. Another option is to 
measure a reduced shear correction rather than calculate it. Equation (fT3l tells us that the needed 
correction is an integral over all modes of the convergence bispectrum, B K . Hence, one could 
conceivably use measurements of shear bispectra to infer the correction and then convert reduced 
shear spectra to shear spectra, thereby recovering the offset-linear scaling form. Bispectra should 
be readily available from weak lensing sur veys since they are useful in their own right for studying 
non-Gaussianity and non-linear clustering (|Cooray & Hull2001bl : lTakada & Jainll2004() . 



5. Conclusions and Caveats 

Reduced shear, the lensing observable measured by averaging galaxy ellipticities, is typically 
approximated by shear for calculational convenience. In this paper, I have explored the error due to 
this approximation in the context of cosmic shear power- and cross- spectra. I have presented the 
general formula, (PT31) . for converting shear spectra to reduced shear spectra, and th e 2D version of 



this p erturbative correction was found to agree with simulations in previous work (|Dodelson et al 



20061) . The correctional term is a projection of the matter bispectrum, which is enhanced by non- 
linear gravitational clustering; the correction is therefore most significant on small scales and for 
highly clustered matter fields, for instance when Og is large. 

I have investigated the extent to which ignoring reduced shear will bias 3D weak lensing 
methods - methods that measure sheared galaxies at multiple redshifts in order to extract informa- 
tion about the expansion of the Universe and the growth of large-scale structure. Weak lensing 
surveys that constrain dark energy parameters via tomography can neglect reduced shear in the 
short term because it is unlikely to be the dominant theoretical systematic error. However, as the 
finer details of shear maps come into focus - observationally and theoretically - reduced shear will 
almost certainly need to be accounted for, as I have demonstrated using a Fisher matrix analysis. 
I have also shown that the shear ratio method known as offset-linear scaling, which probes dark 
energy by isolating the geometric information in cosmic shear, becomes invalid on small angular 
scales because of reduced shear. This method will probably need to be adjusted in order to take 
advantage of information on the smallest angular scales of shear maps. 

The tables in this paper that compare biases in cosmological parameters to their forecasted 
uncertainties are only intended to be a rough guide - they should not be considered definitive 
indicators of the importance of reduced shear, particularly since a real survey will not have the 
perfectly Gaussian likelihood functions assumed by a Fisher matrix analysis. The actual charac- 
teristics of the shear surveys I considered will also determine the importance of the discrepancy; in 
particular, photometric redshift uncertainties will degrade parameter constraints and could change 
the significance of bias. Moreover, combining weak lensing results with different priors or with 
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other dark energy probes (such as type la supernovae) could make a parameter bias better or worse, 
depending on the direction of the bias in parameter space and any tension between the likelihood 
functions being combined. Lastly, the actual difference between shear and reduced shear spectra 
depends strongly on the underlying matter power spectrum, as illustrated by the sensitivity of the 
correction to Cg. 

As a final note, cosmic shear is not the only lensing observable to which reduced shear is 
applicable. Most notably, shear ratio methods are powerful when applied to galaxy-shear cross 
correlations. On small angular scales, we should expect galaxy-shear cross spectra to deviate 
from the offset-linear scaling form in a way similar to the shear- shear deviations plotted in Figure 
HI Furthermore, shear b ispectra and shear- galaxy bispectra may eventually require a 4th order 



reduced shear correction (IDodelson & Zhangll2005l) 



I am grateful to Michael Turner and Scott Dodelson for their guidance on matters both sci- 
entific and professional and to Michael for encouraging me to see "the big picture." Thanks also 
to Bruce Winstein and Jeff Harvey for being on my Ph.D. committee and ensuring that I grad- 
uate with a quality dissertation. I was lucky to have Nobuko McNeill looking out for me at the 
University of Chicago, making my graduate career go as smoothly as possible. Thanks to Neal 
Dalai for suggesting that I take a look at reduced shear in shear ratio methods, and to David Bacon, 
Robert Crittenden, Simon DeDeo, Felipe Marin, Fabian Schmidt, and Martin White for useful dis- 
cussions. Thanks also to Andrew Zentner for the use of his matter power spectrum code, and to 
Jeremy Tinker for his assistance with galaxy bias calculations that ultimately were not included in 
this work. 

This work was supported in part by the Department of Energy and by the Kavli Institute for 
Cosmological Physics through grants NSF PHY-0 114422 and NSF PHY-0551142 and an endow- 
ment from the Kavli Foundation. CS is currently supported by the UK Science and Technology 
Facilities Council. 



-25- 



REFERENCES 

Albrecht, A. et al. 2006, ArXiv Astrophysics e-prints, astro-ph/0609591 

Bacon, D. J., Refregier, A. R., & Ellis, R. S. 2000, MNRAS, 318, 625, arXiv:astro-ph/0003008 
Bartelmann, M., & Schneider, R 2001, Phys. Rep., 340, 291, arXiv:astro-ph/9912508 
Benjamin, J. et al. 2007, MNRAS, 381, 702, arXiv:astro-ph/0703570 
Blandford, R. D., Saust, A. B., Brainerd, T. G., & Villumsen, J. V. 1991, MNRAS, 251, 600 
Cooray, A., & Hu, W. 2001a, ApJ, 554, 56, arXiv:astro-ph/00 12087 

. 2001b, ApJ, 548, 7, arXiv:astro-ph/0004151 

. 2002, ApJ, 574, 19, arXiv:astro-ph/0202411 

Dick, J. 2006, DETFast, http://www.physics.ucdavis.edu/DETFast/ 

Dodelson, S. 2003, Modern cosmology (Modern cosmology / Scott Dodelson. Amsterdam (Nether- 
lands): Academic Press. ISBN 0-12-219141-2, 2003, XIII + 440 p.) 

Dodelson, S., Kolb, E. W., Matarrese, S., Riotto, A., & Zhang, P. 2005, Phys. Rev. D, 72, 103004 

Dodelson, S., Shapiro, C, & White, Martin J., . 2006, Phys. Rev., D73, 023009, astro-ph/0508296 

Dodelson, S., & Zhang, P. 2005, Phys. Rev. D, 72, 083001, arXiv:astro-ph/0501063 

Dunkley, J. et al. 2008, ArXiv e-prints, 803, 0803.0586 

Eisenstein, D. J., & Hu, W. 1999, ApJ, 511, 5, arXiv:astro-ph/97 10252 

Hamana, T., Colombi, S. T., Thion, A., Devriendt, J. E. G. T., Mellier, Y., & Bernardeau, F. 2002, 
MNRAS, 330, 365, arXiv:astro-ph/00 12200 

Heymans, C., & Heavens, A. 2003, ArXiv Astrophysics e-prints, astro-ph/03 10495 

Hirata, C. M., & Seljak, U. 2004, Phys. Rev. D, 70, 063526, arXiv:astro-ph/0406275 

Hu, W. 1999, ApJ, 522, L21, arXiv:astro-ph/9904153 

. 2002, Phys. Rev. D, 66, 083515, arXiv:astro-ph/0208093 

Hu, W., & White, M. 2001, ApJ, 554, 67, arXiv:astro-ph/00 10352 

Huterer, D., & Takada, M. 2005, Astroparticle Physics, 23, 369, arXiv:astro-ph/0412142 



-26- 



Jain, B., & Taylor, A. 2003, Physical Review Letters, 91, 141302, arXiv:astro-ph/0306046 
Kaiser, N. 1992, ApJ, 388, 272 

Kaiser, N., Wilson, G., & Luppino, G. A. 2000, ArXiv Astrophysics e-prints, astro-ph/0003338 

Krause, E. 2008, private communication, paper in preparation 

Ma, Z., Hu, W., & Huterer, D. 2006, ApJ, 636, 21, arXiv:astro-ph/0506614 

National Aeronautics and Space Administration. 2008, JDEM Project Webpage, 
http://universe.nasa.gov/program/probes/jdem.html 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in 
FORTRAN. The art of scientific computing (Cambridge: University Press, — cl992, 2nd 
ed.) 

Rimes, C. D., & Hamilton, A. J. S. 2005, MNRAS, 360, L82, arXiv:astro-ph/0502081 

Rudd, D. H., Zentner, A. R., & Kravtsov, A. V. 2008, ApJ, 672, 19, arXiv:astro-ph/0703741 

Schneider, P., van Waerbeke, L., & Mellier, Y. 2002, A&A, 389, 729, arXiv:astro-ph/0 112441 

Scoccimarro, R., & Couchman, H. M. P. 2001, MNRAS, 325, 1312, arXiv:astro-ph/0009427 

Shapiro, C., & Cooray, A. 2006, Journal of Cosmology and Astro-Particle Physics, 3, 7, 
arXiv : astro-ph/060 1226 

Smith, R. E. et al. 2003, MNRAS, 341, 1311, arXiv:astro-ph/0207664 

Takada, M., & Jain, B. 2004, MNRAS, 348, 897, arXiv:astro-ph/0310125 

Taylor, A. N., Kitching, T. D., Bacon, D. J., & Heavens, A. F. 2007, MNRAS, 374, 1377, 
arXiv : astro-ph/06064 1 6 

Vale, C, & White, M. 2003, ApJ, 592, 699, arXiv:astro-ph/0303555 

Van Waerbeke, L. et al. 2000, A&A, 358, 30, arXiv:astro-ph/0002500 

Wittman, D. M., Tyson, J. A., Kirkman, D., Dell' Antonio, I., & Bernstein, G. 2000, Nature, 405, 
143, arXiv:astro-ph/0003014 

Zhang, J., Hui, L., & Stebbins, A. 2005, ApJ, 635, 806, arXiv:astro-ph/03 12348 



This preprint was prepared with the AAS I^TgX macros v5.2. 



